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PERTURBATIVE QCD- AND POWER-CORRECTED HADRON SPECTRA AND SPECTRAL 
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A. ALI, G. HILLER 
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Leading order (in a s ) perturbative QCD and power (1/m^) corrections to the hadronic invariant mass and hadron energy spectra in 
the decay B — » X s £ + £~ are reviewed in the standard model using the heavy quark expansion technique (HQET). In particular, the 
first two hadronic moments (Sjj) and (E^), n = 1,2, are presented working out their sensitivity on the HQET parameters Ai and 
A. Data from the forthcoming B facilities can be used to measure the short-distance contribution in B — » X s £ + £~ and determine 
the HQET parameters from the moments {S^). This could be combined with complementary constraints from the decay B — > X£i>£ 
to determine these parameters precisely. 



1 Introduction 

The semileptonic inclusive decays B — > X£ + £~ , where 
£ ± = e ± ,/i ± ,r ± and X represents a system of light 
hadronic states, offer, together with the radiative electro- 
magnetic penguin decays B — > A + 7, presently the most 
popular testing grounds for the standard model (SM) in 
the flavour sector. In this contribution, we summarize 
the main steps in the derivation of the hadron spectra 
and hadron spectral moments in B —>■ X s £ + £~ using per- 
turbative CDCD and the heavy quark expansion technique 
HQET Sffl, published recently by us ii. This work, 
which incorporates the leading order (in a s ) perturbative 
QCD and power (l/m|) corrections to the hadronic spec- 
tra, complements the derivation of the dilepton invariant 
mass spectrum ancLthe forward-backward asymmetry of 
the charged leptontl, calculated in the HQET framework 
some time ago by us in collaboration with T. Morozumi 
and L. Handokolil. (See, also Buchalla ct al. 13.) Both the 
hadron and dilepton spectra are needed to distinguish the 
signal (B — > X s £ + £~ ) from the background processes and 
in estimating the effects of the experimental selection cri- 
terion. We shall concentrate here on the short-distance 
contribution which can be extracted from data with the 
help of judicious cuts, such as those employed recently by 
the CLEO collaborationtl The residual effects from the 
resonant (long-distance) contributions have been studied 
in these distributions elsewhere EHMJ, to which we refer 
for details and references to the earlier work. 

We also underline the theoretical interest in measur- 
ing the first few hadronic spectral moments (Sjj) and 
(E%) (n = 1,2). The former are sensitive to the HQET 
parameters A and Ai; we work out this dependence nu- 
merically and argue that a combined analysis of the mo- 
ments and spectra in B — > X 8 £ + £~andi B — > X£vt will 
allow to determine the HQET parameters with a high 
precision. Since these parameters are endemic to a large 
class of phenomena in B decays, their precise knowledge 



is of great advantage in reducing the theoretical errors in 
the determination of the CKM matrix elements Vtdi Vts, 
V cb and V u b- 

2 Kinematics and HQET Relations 

We start with the definition of the kinematics of the de- 
cay at the parton level, b(pb) — > s(p s )(+g(p g ))+£ + (p + ) + 
£~(p-), where g denotes a gluon from the 0(a s ) correc- 
tion. The corresponding kinematics at the hadron level 
can be written as: B(pb) — > X s (Ph) +^ + (p+) + 
We define by q the momentum transfer to the lepton pair 
q = p + + p- and s = q 2 is the invariant dilepton mass 
squared. We shall also need the variable u defined as 
u = —(pb — P+) 2 + (pb — P-) 2 ■ The hadronic invariant 
mass and the hadron energy in the final state is denoted 
by Sh and Eh, respectively; corresponding quantities 
at parton level are the invariant mass sq and the scaled 
parton energy xq = From energy- momentum con- 
servation, the following equalities hold in the 6-quark, 
equivalently i?-meson, rest frame (y = (1,0,0,0)): 

Xq = 1 — V ■ q , s = 1 — 2v ■ q + s , 
E H =m B —v-q, Sh = m B - 2m B v ■ q + s , (1) 

where dimensionless variables with a hat are scaled by 
the 6-quark mass, e.g., s = -^r, rh s = — etc. Here, the 
4-vector v denotes the velocity of both the &-quark and 
the £?-meson, pb = mbV and ps = m B v. 

The relation between the i?-meson and 6-quark mass 
is given by the HQET mass relation m B = nib + A — 
l/2r7ib(Ai + 3A2) + . . ., where the ellipses denote terms 
higher order in l/mb. The quantity A2 is known precisely 
from the B* — B mass difference, with A2 — 0.12 GeV 2 . 
The other two parameters are considerably uncertain at 
present and are of interest here. 

The hadronic variables E B and Sh can be expressed 
in terms of the partonic variables xq and §0 by the 
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following relations 



with 



T Ai + 3A 2 / j Ai + 3A 2 
E H = A L - + [m B - A + -i- I x 



2m b \ 2m B 

Sff = m 2 s + A 2 + (m| - 2 Am 5 + A 2 + A : + 3A 2 ) (s - ™s) 
+ (2Am B - 2A 2 - Ai - 3A 2 )a;o + . . . . 

The dominant non-perturbative effect on the hadron 
spectra is essentially determined by the binding energy 
A = tub —mb + in terms of which one has the following 
transformation: 



E -> E H = A + E + . . . , 
so -> S H = s + SA^q + A 2 



(2) 



Thus, changing the variables of integration (sq,Eq) — > 
(sh, E ) and integrating over E in the range s/Sh — A < 
Eq < 1/2itlb{Sh — 2Am B + m B ), one gets an invariant 
hadron mass spectrum dT/dSn in the kinematic range 
A 2 < Sh < m B . In particular, already for the partonic 
decay b — > s£ + £~ with m s = 0, and hence so = 0, one gets 
a non-trivial distribution in Sh for A 2 < Sh < Arris- 
The kinematic boundary of the distribution dT/dSn is 
extended by the bremsstrahlung process b — > s + g + 
£ + £~, where now Arris < Sh < rn B (with m s = 0). 
The 0(a s ) contribution leads to a double logarithmic 
(but integrable) singularity at Sh — Am b ■ Perturbation 
theory is valid for A 2 < Sh < m \i A 2 > Arris- 

3 Matrix Element for B -> X s £+£- in the Effec- 
tive Hamiltonian Approach 

The effective Hamiltonian governing the decay B — > 
X s £ + £~is given aslil: 



4gj 
"V2 



e 



(3) 



V Ci(/i)Oi + C 7 {ri)-— i s a a^{m b R + m s L)b a F^ 

2 — 1 

e 2 - e 2 - 

lbn z loir z 



where Gf is the Fermi coupling constant, L(R) = 1/2(1=F 
75), and C; are the Wilson coefficients. Note that the 
chromo-magnetic operator does not contribute to the de- 
cay B — ► X s £ + £~in the approximation which we use here. 

The matrix element for the decay B — > X s ^ + ^~can 
be factorized into a leptonic and a hadronic part as 



M(B^xj+e-) = ^v*v tb (r L „L L " 

V 2n 



T R ,L Rfi 



(4) 



L L ' R ^ =£^L(R)£ 

V L/R = 



RlJcf(s) T C w + 2cf^ 



+2m s C 7 cff 7 M |i 



b. 



(5) 



(6) 



with Cf ff = C 7 -C 5 /3-C 6 - The effective Wilson coef- 
ficient C| ff (s) receives contributions from various pieces. 
The resonant cc states also contribute to Cg (s); hence 
the contribution given below is just the perturbative part: 



Cf ff (S)=C 9 » 7 (M + r(.s). 
i) repres 



(7) 



Here 77 (s) and Y(s) represent, respectively, the 0(a s ) 
correction HI and th e on^ loop matrix element of the 
Four-Fermi operators!; 

With the help of the above expressions, the differen- 
tial decay width becomes on using p± = (E± , p± ) , 



dr 



1 G F 2 a 2 
2m r 2tt 2 



KVa, 



2 d 3 p + d 3 p_ 
(2n) 3 2E+ (2n) 3 2E„ 



(8) 



where W^: R and L^;, fl are the hadronic and leptonic tcn- 
sors, respectively, and can be seen in the literature III. 

L / R. 

The hadronic tensor W^l is related to the discontinuity 

L I R 

in the forward scattering amplitude, denoted by T^l , 
through the relation W^J R — 2lmT^J R . Transforming 
the integration variables to s, u and v ■ q, one can express 
the triple differential distribution in B — > X s £ + £r&s: 



dr 



1 Gp 2 a 2 mj 



dudsd(vg) 2m B 2vr 2 256tt 4 ' V * s ^* b ' 



x 21m 



(9) 



Using Lorentz decomposition, the tensor T^,, can be ex- 
panded in terms of three structure functions 7$, 



rpL/R 



-T, 



L/R 



9^ 



T, 



L/R 



Vn v v + T^ /R ie^afj v a if , 



(10) 

where the ones which do not contribute to the amplitude 
for massless leptons have been neglected. 

4 Hadron Spectra in B -> X s £+£~ 

We discuss first the perturbative 0(a s ) corrections re- 
calling that only the matrix element of the operator 
Og = e 2 /(167r 2 )s Q 7 M L6 Q €7 M € is subject to such correc- 
tions. The corrected hadron energy spectrum in B — > 
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X s £ + £ can be obtained by using the existing results in 
the literature on the decay B — > Xlvi by decomposing 
the vector current in 9 as V = (V - A) /2 + (V + A) /2. 
The (V — A) and fV + vl) currents yield the same hadron 
energy spectrum 113 and there is no interference term 
present in this spectrum for massless leptons. So, the 
correction for the vector current case can be taken from 
the corresponding result for the charged (V — A) caseEL?J. 

The 0(a s ) perturbative QCD correction for the 
hadronic invariant mass is discussed next. As already 
mentioned, the decay b — > s + £ + + l~ yields a delta 
function at so = m 2 and hence only the bremsstrahhmg 
diagrams b — > s+g+£ + +£~ contribute in the range m 2 < 
Sq < 1. The resulting distribution dB(B — > X s £ + £~)/dso 
in the parton model in the 0(a s ) approximation and 
thejSudakov exponentiated form can be seen in our pa- 
per B. We remark that the Sudakov exponentiated dou- 
ble differential distribution for the decay B — > X u £vt 
has been derived by Greub and ReyO, which we have 
checked and used after changing the normalization for 
B — > X s £ + £~. The hadronic invariant mass spectrum 
dB{B -► X s £+£-)/dS H , shown in_Fig. |l| depends rather 
sensitively on mt (or equivalently A). An analogous anal- 
ysis for the decay B — > X u £vg has beeiLperformed earlier, 
with very similar qualitative results LL2I. 

Next, we discuss the power corrections to the 
hadronic spectra. The structure functions T^ R in the 
hadronic tensor in Eq. ( |Io| ) can be expanded in inverse 
powers of with the help of the HQET techniquesllrH'B. 
The leading term in this expansion, i.e., 0(m®), repro- 
duces the parton model result HD'HU. In HQET, the next 
to leading power corrections are parameterized in terms 
of Ai and A2. After contracting the hadronic and lep- 
tonic tensors and with the help of the kinematic identi- 
ties given in Eq. (Q) , we can make the dependence on xq 
and so explicit, 

L L "*" = m b 2 {2(1 - 2xo + so)T x L l R 

T 2 L ^T(l-2x + s )uT 3 L / R Y 

(11) 



-u - s 



By integrating Eq. (^) over u, the double differential 
power corrected spectrum can be expressed asll: 



d 2 B 
dxo dso 



8 I 

-B a lm\ x 2 - s {(1 - 2a; + s )T 1 (s ,x ) 

7T v 



X - S r 



T2(s ,x ) > + 0(X l a s ) . 



(12) 



The structure function T3 does not contribute to the dou- 
ble differential distribution and we do not consider it any 
further. The functions Ti(sq, xq) and T2(sq, xo), together 



m b = 4.8 GeV 
m b =4.7 GeV 
m h =4.9 GeV 




3 4 
S H [GeV ! ] 



Figure 1: The differential branching ratio dB(B -* X s £+£-)/dS H 
in the hadronic invariant mass, Sh, shown for three values of m;, 
in the range where only bremsstrahhmg diagrams contribute. 



with other details of the calculations, have been given by 
us elsewhere H. 

The branching ratio for B — > X s £ + £~is usually ex- 
pressed in terms of the measured semileptonic branching 
ratio B s i for the decay B — > XJLvi. This fixes the nor- 
malization constant Bq to be, 



Bo = B s 



3a 2 \V£V tb \' 



1 



16^ 2 \V cb \ 2 f(m c ) K (m c ) ' 



(13) 



where f(rh c ) is the phase space factor for T(B — > X c £vi) 
and «(m c ) accounts for both the 0(a s ) QCD correction 
to the semileptonic decay widthllE and the leading order 
(1/mfc) 2 power correctionlll. The hadron energy spectrum 
can now be obtained by integrating over so with the kine- 
matic boundaries: max(rh 2 , —1 + 2xq + Amf) < sq < Xq, 
iris < a^o < 5(1 + ™s — 4m 2 ). The hadron energy spec- 
trum dB(B — > X s £ + £~)/dEo in the parton model (dotted 
line) and including leading power corrections (solid line) 
are shown in Fig. ^. For rrib/2 < Eq < m& the two distri- 
butions coincide. Note that the l/m 2 -expansion breaks 
down near the lower end-point of the hadron energy spec- 
trum and at the cc threshold. Hence, only suitably aver- 
aged spectra are useful for comparison with experiments 
in these regions. Apart from these regions, the HQET 
and parton model spectra are remarkably close to each 
other. 

5 Hadron Spectral Moments in B — > X s £ + £~ 

The lowest spectral moments in the decay B — > 
X s l + £~a,t the parton level are worked out by taking into 
account the two types of corrections discussed earlier, 
namely the leading power l/m,b and the perturbative 
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Figure 2: Hadron energy spectrum dB(B — ► X s iP~l~)/dEo in the 
parton model (dotted line) and including leading power corrections 
(solid line). For m.(,/2 < Eq < ra^ the two distributions coincide. 



0{a s ) corrections. To that end, we define the moments 
for integers n and m: 



M 



(n,m) 

l+l- = 



1 



AB 



(s -m 2 s rx^—— 
B J ds a dx 



dsndiro , (14) 



which obey (x™(§o — mj) n ) = ^-.M[™j_ . These mo- 
ments can be expanded as a double Taylor series in a s 
and 1/mb- 



(n,m) 



(n,m) • _ 



(15) 
0, 1,2, into 



with a further decomposition of D 
pieces from different Wilson coefficients: 

D (n,m) = a (n,m) C eS 2 + ^ + jM C f + . 

(16) 

The terms -y| n,m ' and S^ n in Eq. (|l6| ) result from 
the terms proportional to i?e(Cf ff )Cf ff and |Cf ff | 2 in 



Eq 

J 



respectively. The explicit expressions for 



AtL,m) (n,m) f (n,m) 

p. , 7> , o - are given m our paper L 



The leading perturbative contributions for the 
hadronic invariant mass and hadron energy moments can 
be obtained analytically, 

A (o,o) = 25 ~ 47f2 A (i,o) = a(2,0) = J_ 

9 675 ' 486 ' 

A (o,i) _ 1381 - 2107T 2 ^ (0|2) _ 2257- 320^ 2 ^ 



1350 



5400 



The zeroth moment n — m — is needed for the normal- 
ization; the result for A^ ' ^ was first derived by Cabibbo 



and MaianillJI. Likewise, the first mixed moment A^ 1 ' 1 ) 
can be extracted from the results for the decay B — > Xlvi 
c3 after changing the normalization, A^ 1,1 ) = 3/50. For 
the lowest order parton model contribution Dq ' m \ we 
find, in agreement withE^, that the first two hadronic 



rn 

■„2 



), ((so — m 2 ) 2 ) and the 
)) vanish: D^ 0) = 0, for 



invariant mass moments (§o - 
first mixed moment (xq(sq — 
n = 1,2 and D^' X) = . 

Using the expressions for the HQET moments de- 
rived by us El, we present the numerical results for the 
hadronic moments in B — ► X s £ + £~. The parame- 
ters used are : m s = 0.2 GeV,TO c = 1.4 GeV, nib = 
4.8 GeV,m t = 175 ± 5 GeV,/x = m fe +™\ 2 , a s (m z ) = 
0.117 + 0.005 , a" 1 = 129. We find for the short-distance 
hadronic moments, valid up to 0(a s /m 2 B , l/m B ): 



(Sh) =m 2 B (^- + 0.093 



0.069 



A a s 



TOfl 7T 



A A 2 A A 

+0.735 + 0.243— + 0.273^- - 0.513— §-) . 



m B 



(S 2 H ) =mt(0.0071— 

TT 



0.138 

m B tt 



)>> r 



+0.587- 



A 2 



(Eh) = 



+0.012- 



A 2 



0.196 



0.024 



(E H ) 



0.147m|(l + 0.324 
A 2 



+2.740- 



0.299 



■4) 

m B 




(18) 




0-352-5-2- . 


1- 1.691 — 


TT 


m B n 


m B 




+ 1.07(A), 




m B 


m% 




a s 


0.128 A*. 


V 2.954— 


n 


tob tt 






+ 0.162 A), 




m B 


m B 





where the numbers shown correspond to the central val- 
ues of the parameters. 

The dependence of the hadronic moments given in 
Eq. ( p"8| ) on the HQET parameters Ai and A has been 
worked out numerically. In doing this, the theoretical 
errors on these moments following from the errors on the 
input parameters mt, a s and the scale [i have been esti- 
mated by varying these parameters in the indicated ±lcr 
ranges, one at a time, and adding the individual errors 
in quadrature. The correlations on the HQET param- 
eters Ai and A which follow from (assumed) fixed val- 
ues of the hadronic invariant mass moments (Sh) and 
(S 2 H ) (calculated using A = 0.39 GeV, Ai = -0.2 GeV 2 
and A 2 = 0.12 GeV 2 ) are shown in Fig. | (for the de- 
cay B — ► A" s /i + /.t~). The (Ai-A) correlation from the 
analysis of Gremm et al. for the electron energy spec- 
trum in B — > X£vj> is shown as an ellipse in this fig- 
ure. With the measurements of (Sh) and (S H ) in the 
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Figure 3: (Sff) (solid bands) and (S^) (dashed bands) correlation 
in (Ai-A) space for the decay B — > X s f+£ — . The correlation from 
the analysis of the decay B — > X£u£ by Gremm et al. lij is shown 
as an ellipse. 



decay B — ► X s £ + £~, one has to solve the experimen- 
tal numbers on the l.h.s. of Eq. ( |l8| ) for A x and A. It 
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in B — > X s £ + £~ remain to be defined, but for the cuts 
used by the CLEO collaboration we have studiedJiie ef- 
fects in the HQET-like Fermi motion (FM) modelcil. We 
find that the hadronic moments in the HQET and FM 
model are very similar and CLEO-tyne cuts remove the 
bulk of the cc resonant contributions B. 

In summary, we have calculated the dominant con- 
tributions to the hadron spectra and spectral moments 
in B — ► X s £ + £~ including contributions up to terms of 
0{a s /m 2 B , l/m%). We have presented the results on the 
spectral hadronic moments (Eg) and (S H ) for n = 1,2 
and have worked out their dependence on the HQET pa- 
rameters A and Ai. The correlations in B — > X s £ + £~ are 
shown to be different than the ones in the semileptonic 
decay B — > X£v^. This complementarity allows, in prin- 
ciple, a powerful method to determine them precisely 
from data on B — > X£vi and B — > X s £ + £~m forthcoming 
high luminosity B facilities. 
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